Mon. Not. R. Astron. Soc. 000, [l]{T5] (2010) Printed 8 February 2012 (MN stylo file v2.2) 



A hybrid steady-state magnet ohydro dynamic dust-driven 
stellar wind model for AGB stars 



O 



6 



Anand Thirumalai^* and Jeremy S. Heyl^* 

^ University of British Columbia, 6224 Agricultural Road, Vancovuer, British Columbia, V6T IZl, Canada 



8 February 2012 



ABSTRACT 

We present calculations for a magnetised hybrid wind model for Asymptotic Giant 
Branch (AGB) stars. The model incorporates a canonical Weber-Davis (WD) stellar 
wind with dust grains in the envelope of an AGB star. The resulting hybrid picture 
preserves traits of both types of winds. It is seen that this combination requires that 
the dust-parameter {Td) be less than unity in order to achieve an outflow. The emer- 
gence of critical points in the wind changes the nature of the dust-driven outflow, 
simultaneously, the presence of a dust condensation radius changes the morphology 
of the magnetohydrodynamic (MHD) solutions for the wind. In this context, we ad- 
ditionally investigate the effect of having magnetic-cold spots on the equator of an 
AGB star and its implications for dust formation; which are seen to be consistent 
with previous findings. 

Key words: AGB wind, Weber-Davis wind, MHD, dust-driven wind. 
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1 INTRODUCTION 

In recent years several observations of evolved stars and 
planetary nebulae have indicated that magnetic fields may 
exist in these objects (e.g. Amiri et al.||2010l [Jordan e t al.| 
[2005 , Herpin et al. 2006 [ [ Sabin et al.[20 07a Sabin ct al. 2007V 
and [Miranda et al.||2001 l. The inferred magnetic fields in 
such objects, at the distances of the masers, indicate a va- 
riety of field strengths ranging from a few milligauss to a 
hundred gauss or so ( jVlcmmings et al.|2005[[2C)06, ). This in 
turn, may indicate a variety of field strengths at the sur- 
faces of such objects. The idea of studying the effects of a 
magnetic field on the nature of the AGB wind has been car- 
ried out mainly with regard to explaining the diversity of 



the observed shapes (see for example Chu et al. 1987 and 
[Stanghellini et al.|1993[ and references therein) of planetary 
nebulae. As a result, the modelling of magnetic fields in these 
stars has focussed on the final stages of the AGB phase, at 
the very tip of the AGB, or on the post-AGB phase itself, 



employing the so called interacting wind scenario (see Frank 
pM9l[Franco et al.|200ll|Gardiner fc Frank|2001[ for outlines 
of the key points of interacting wind models). 

It has therefore been argued that magnetic fields play 
a vital role in shaping planetary nebulae by several re- 
searchers. However, it is not completely certain whether 
magnetic fields play a dynamic global or local role in this re- 
gard. The reader is referred to work carried out by Sober and 
co-workers (e.g. [Soker fc Kastner|20"03| [Soker fc Zoabi.2002[ 
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[Soker fc Cl8yton|[T999l [Soker fc Harpai[[T999| 

for discussion regarding this topic. Moreover, 
current state-of-the-art MHD models for the winds of stars 
at the tip of the AGB, do not take into account the effect of 
radiation pressure on the dust grains in the envelope of an 
AGB star. It is generally thought that the mass loss in AGB 
stars is largely governed by this mechanism coupled with 



strong stellar pulsations ( 


Lamers fc Cassinelli|1999, ^Elitzur 


et al.[[2003[ [Bowen[[1988[ 


Bowen & Willson[[1991[ jFleischer 


et al.[[1992[ Wood|1979 Bedijn[[1988 


\. However, to the best 



of our knowledge, there have not been any investigations 
combining a standard dust-driven wind scenario with MHD 
effects, in the literature. Such a study, given the importance 
of magnetic activity in AGB stars, would help bring together 
two different sub-classes of stellar winds. This is the aim of 
the current paper; to investigate the implications of com- 
bining magneto-rotational effects with a dust-driven wind 
in AGB stars. It is to be mentioned that this would be ap- 
plicable at the early stages of the AGB phase, long before 
the interacting wind scenario becomes important, wherein 
the models mentioned above would be more likely candi- 
dates for describing the outflows. 

At this juncture, we conduct a brief survey of the avail- 
able literature with regard to both dust-driven winds and 
magneto-rotational equatorial winds. However, the reader 
is referred to a review by Tsinganos (2007 and refer- 



ences therein) for MHD outflows and likewise, reviews by 



Dorschner[ | [2003[ ) and [Habing fc 01ofsson| ( [2004[ ) for histori- 
cal reviews of AGB stars and dust-driven winds. Thereafter, 
the model developed in the current study, is elucidated. 
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With regard to magneto-rotational equatorial winds, 



seminal work was carried out by Weber & Davis ( 1967 
WD). They formulated a steady state description for the 
radial and azimuthal components of the solar wind's mo- 
mentum. Essentially the same results were also arrived at 



independently, by Mestel (19671. These studies formed the 



groundwork upon which further studies were conducted. 
Thereafter, magnetic braking by a stellar wind was also in- 



vestigated by Mestel ( 1968 1 and later in greater detail by 
Okamoto ( 1974 1975 1, wherein the theory was extended to 
cover a variety of field configurations with poloidal fields. 
The Weber and Davis equatorial wind theory was extended 



by Goldreich & Julian ( 1970 1 towards a relativistic treat- 



ment of the wind, including the effects of pressure and grav- 
itation. Michel ( 1969 1 carried out a similar analysis of the 
WD model but neglected pressure and gravity and thus, rel- 
ativistic magnetosonic critical points do not appear in his 
model. The first effort to investigate the importance of the 
interactions between the gas and the magnetic field, with 
regard to determining properties of the structure and dy- 
namics of the solar corona, was carried out by |Pneuman fc| 
Kopp ( 1971 1. They found that the assumed dipolar field con- 



figuration had a profound effect on the solar wind in creating 
streamers. |Yeh| ( |1976[ ) conducted a parametric study of the 
mass and angular momentum effluxes of magneto-rotational 
stellar winds and found that the mass efflux would be large, 
if the mass of the star was small, with a large radius, pro- 
vided the stellar corona was dense and hot. Simultaneously, 
he found that the angular momentum efflux became greater 
when the magnetic field and stellar rotation parameters were 



increased. Belcher & MacGregor ( 1976 1 revisited the Weber 



and Davis theory and identified two regimes; the so-called 
slow and fast magnetic rotators, which are defined by the 
ratios of the Michel and Parker velocities of the wind. The 
former is related to the ratio of the magneto-rotational flux 
to the mass flux and the latter is related to the squares of 
the sound speed, escape velocity and radial bulk gas veloc- 
ity at the surface. Presently, the reader is referred to Belcher^ 
|fc MacGregor] ( |1976[ ) for mathematical expressions for these 
quantities. Belcher and MacGregor were able to delineate 
the angular momentum evolution of the two types of rota- 
tors on the main sequence. [Barker fc Marlborough] ( [19821 ) 
extended the Weber-Davis theory for the case of non-zero 
photospheric mass loss, showing that including a dimen- 
sionless constant corrected the original theory. The original 
WD theory was also re-analysed elsewhere, with an effort to 
study standing MHD shocks by 'Chakrabarti (19901. Vari- 



ous solutions of the WD solution topology were studied and 
he found that many of them allowed for MHD shock forma- 
tion, in both accretion and winds, of a compact magnetised 
object. With the advent of increased computational capabil- 
ities of computers in the 1980's, Sakurai ( 1985') successfully 



generalised the original Weber and Davis theory to two di- 
mensions. There appear the usual slow and fast modes in his 
solution, as in the original Weber and Davis theory; however, 
in this case, the momentum equation was found to be singu- 
lar on an Alfven surface and regularising the solution on this 
surface alongside the boundary condition at the photosphere 
uniquely determined the solution of the two-dimensional 
problem. Keppens & Goedbloed ( 1999 2000 1 have over the 



past few years, developed two- and three-dimensional MHD 
models for investigating magento-rotational stellar winds. 



Their models have clearly shown that the dipolar nature of 
the magnetic fleld structure is important for stellar winds. 
At the same time they have shown that the poloidal com- 
ponent leads to density enhancement along the equatorial 
region. They also flnd that their trans-sonic wind solutions 
have dead zones which have a latitudinal dependence; this 
can be traced back to the configuration of the magnetic field. 
They have also investigated shock formation in magneto- 
rotational outflows ( Zaliznyak et al.||2003 see) and the na- 
ture and formation of Kelvin-Helmholtz instabilities (e.g. 
Keppens et al.1'1999'). These studies indicate that the mass 



efflux from magneto-rotational outflows can be asymmet- 
ric in nature. Their models also have been used to study 
the evolution of rotational velocity distribution in late-type 
stars while on the main sequence. The stellar winds in these 
stars were assumed to be WD-like steady-state winds (see 
Keppens et al.|1995 l. 



Concomitant to the development of MHD winds, the 
field of dust-driven winds in evolved stars has flourished as 
a separate research focus altogether. The early models of 
dust-driven outflows from AGB stars were one-dimensional 
and did not include radiative transfer through the dust- 



laden envelope and dynamical dust formation (see Lamers 
fc Cassinelli[l999| for a brief review). Later one-dimensional 



models involved dynamics of dust formation and growth in 
a time-dependent manner building upon initial studies of C- 



type AGB stars (e.g. 


Gail & Sedlmayr||1988 


|1990| [Fleischer et al. 


|1991|). These models 



due to the dynamical interaction between dust and gas and 
the radiation field, dust-laden shells were formed in the en- 
velopes of AGB stars (e.g |Fleischer et al.|1992| and dynam- 
ical instabilities in the flow also manifested themselves lead- 
ing to aspherical mass loss (se e [Fleischer et al.|1995, Hoefnerl 
|et al.llTgtsl [Dorfl fc Hoefner|[l996| . Models have also been 
constructed that investigate dust grain drift through the gas 
which modifies dust growth rates and the efficiency of the 
wind acceleration process (see Sandin &H6fner 2004'). There 
have also been further improvements in treating the radia- 
tive transfer in a frequency dependent way alongside detailed 
micro-physics of molecular and grain opacities (e.g. |Hofner| 



et al. 2003 Andersen et al. 2003). The reader is presently 



referred to a recent review by Hofner (2008) on the status 



of modern-day radiation-hydrodynamical modelling of dust- 
driven winds of AGB stars. The modelling of dust-driven 
winds has elsewhere, been extended to more than one dimen- 
sion, capturing the aspherical nature of the outflow in a very 
clear way (see|Woitke fc Niccolini '" 2005 , Woitkc 2005] |2006 



Woitke & Quirrenbach 2008 Woitke 2008). These models 



incorporate hydrodynamics with radiation pressure on the 
dust, equilibrium chemistry for the nucleation of dust in 
a time-dependent way, and they take into account radiative 
transfer in a frequency dependent manner. They were able to 
capture in their simulations, many of the highly dynamical 
aspects of the outflow, including turbulent and inhomoge- 
neous dust formation and Rayleigh- Taylor flow instabilities 
that result in cloud-like structure formations in the efflux. 
With the help of such models it has become clear that the 
winds in these stars are far more complicated than simple 
one-dimensional (spherically symmetric) pictures and may 
have an impact on shaping planetary nebulae due to their 
asphericity (e.g. Reimers et al.^2000) or by having an impact 
on the superwind at the end of the AGB, as is suggested by 
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Lagadec & Zijlstra' ('2008') , with regard to the abundance of 2 THE HYBRID WIND MODEL 



carbon in the envelopes of some AGB stars. The latter sug- 
gests that the effort to include complicated micro-chemistry 
of the dust grains, as is done in modern AGB wind models, 
may prove to be a key factor in resolving issues regarding the 
superwind in the late stages of the AGB phase. For a more 
thorough description of AGB envelopes and the stochastic 
nature of dust formation process itself, the reader is referred 
to [Habing fc 01ofsson| (|2004[) and to [Dirks et al.] (|2008[) as 



well as references therein. 



It has also been argued that asymmetric mass loss on 
the AGB may result in kicks and spin being delivered to the 
nascent white dwarf within (see Spruit|1998 1 . More recently, 
the effect of kicks to white dwarfs has been investigated with 
regard to their impact on the dynamics of globular clusters 
and binaries (see [Davis et al.|[2008l [Heyl[[2007aP [2008a|b 



[Heyl fc Penrice|2009[ ). These studies indicate that the study 
of AGB winds in relation to asphericity is rather important. 
In the literature, hybrid equatorial winds have only been 
considered for rotating hot stars (see MacGregor et al.|1992 



MacGregor fc Friend|[T987l [Friend fc Mac Grcgor 1984| [Poe[ 



et al.|198 9 : Poc & Friend 1986 ). MHD stellar winds for AGB 



stars have been simulated for a coupled disk and star system, 
clearly showing that magnetocentrifugal winds can occur in 
planetary nebulae and AGB stars, at the tip of the AGB 
(see Blackman et al. 2001 and references therein). Other 



researchers have considered a hybrid wind for AGB stars by 



coupling Alfven waves with radiation pressure (see Falceta- 



Gongalves fc Jatenco-Pereira|2002} [Falceta-Gongalves et al 



2006 1, showing that it is possible to obtain low- velocity mass 



efflux in supergiant cool stars through such a mechanism. A 
similar study has also been carried out for Wolf-Rayet stars 
(e.g Dos Santos et al.|1993 1, of course obtaining much faster 
winds. From the brief survey of the literature conducted 
above, it is evident that winds from AGB stars whether 
during the initial or final phases of the AGB can have asym- 
metries that result in aspherical mass loss either due to 
purely MHD effects, or due to the dynamics between the 
dust and gas and the radiation field. However, present day 
dust-driven models still do not include magneto-rotational 
effects, particularly when more and more observations indi- 
cate the presence of magnetic fields in these objects and the 
pure MHD wind models do not include the effect of dust con- 
densation and radiation pressure in the envelopes of these 
stars. The aim of the current work is to investigate the impli- 
cations of including magneto-rotational effects with a sim- 
plified dust-driven model in 1.5 dimensions; the azimuthal 
terms are determined entirely from their dependence on the 
radial terms (the standard WD-picture). We shall however, 
for the sake of simplicity, be considering a steady-state case. 
It is to be re-iterated that the current model may only be 
valid at the early phases of the AGB, long before the onset 
of the superwind at the end of the AGB. 



We present the underlying assumptions of the model in 
the following section (!|2|. In ijsjwe present the numerical de- 
tails, thereafter in sQthe results are presented and discussed. 
Finally, in !j5| the paper is summarised and future avenues 
for development of the current model are briefiy discussed. 



We begin with the standard picture of a WD-equatorial 
wind; we shall be following their definitions and general 
derivation closely. However, we shall be including the dust as 
a second fluid. As in the WD model we shall assume com- 
plete axial symmetry and the customary explicit form for 
the magnetic field within the equatorial plane, 



B = Br{r)f + B4r) 



(1) 



while the velocity of the gas and the dust can respectively, 
be written as, 



u = Ur{r)f + u^{r) 
and 

V = Vr(r)f + Vq,{r)i 



(2) 



(3) 



As can be seen, the velocity fields are functions of the radial 
distance alone. In addition, there is no time dependence, 
explicit or implicit, thus conforming with the steady-state 
assumption. Accordingly, the continuity equation for the gas 
and dust combined can be written as. 



pur -I- ndmdvr Q{r — rd) ~ constant. 



(4) 



where, p is the gas density in cgs units, Ud is the number 
density of the dust grains and rud is the mass of a single dust 
grain. Q{r — rd) is the standard Heaviside function which is 
equal to unity for r ^ rd- Prior to the dust formation radius, 
there is only one fiuid, namely the gas. Please note that we 
have dropped the subscript r from the radial velocities for 
brevity and future convenience. We are explicitly assuming, 
for the sake of simplicity, that all dust formation occurs at 
a certain distance (rd) from the centre of the star. Beyond 
this, there is no further condensation of dust. We shall also 
assume that all the dust grains are identical and perfectly 
spherical, with a radi us of a = 0.05/xm and a ma ss density of 
Pd ~ 2.25g/cm^ (e.g. Lamers fc Cassinelli|l999 1. This yields 
the mass for an individual dust grain; rrid = ^na^pd- 

In addition, we shall assume that the dust-to-gas ratio 
in the stellar wind is given by, 



ridrrid 



200 



(5) 



following Lamers & Cassinelli (19991, in other words, 



ndfnd <^ p- The maximum value for this ratio is at the dust 
formation radius and it decreases monotonically thereafter. 
For deriving the momentum equations for the gas and the 
dust, our starting point is the Euler equation, one for each 
of the fiuids. For the gas we can write. 



{u- \7)u+ -Vp + 



GM, 



1 

pc 



J X B 



rid 



/o =0, 



(6) 



where p is the gas pressure and J is the current density. It 
is assumed here implicitly that there is no relative motion 
of the ions with respect to the neutrals. The first and sec- 
ond terms are related to the velocity and pressure gradients 
respectively, the third term is the gravitational acceleration 
on the gas, the fourth term is the Lorentz force divided by 
the gas density and the final term is proportional to the 
drag force that is experienced by the gas, due to the dust 
grains moving through it. Assuming force-free MHD it can 
be shown that (e.g. Weber fc Davis|1967 l 
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r{uB^ — u^Br) = constant = —RgQBrfi, 



(7) 



where, -Ro is the steUar radius (radius of the photosphere), 
Q is the rotation rate of the star and _Br,o is the radial com- 
ponent of the magnetic field at the stellar surface. Requiring 
that V • -B = 0, yields the familiar relation 



r^Br = -Ro^r.O- 



(8) 



Presently, let us turn our attention to the Euler equation 
for the second fluid, viz., the dust grains. The rationale is 
that obtaining expressions for the drag force components 
will allow us to re- write Eq. the Euler equation for the 
gas. The Euler equation for the dust grains can be written 
as, 



,^ GNU , 
{v ■ \)v H r - 



r H 

rud 



fD=0. 



(9) 



In the above equation, Qrp is the radiation pressure mean 
efficiency (e.g. [Lamers fc Cassinelli| [19991 ). current 
study, since our aim is to present a simplistic picture we 
shall not be calculating this term. L« is the luminosity of 
the star and c is the speed of light. The first term in Eq. (|9| 
is related to the velocity gradient of the dust grains, the 
second term is the gravitational acceleration experienced by 
the dust grain, the third term is the radiation pressure that 
the dust grain experiences, that drives it outward and the 
final term is the drag force per unit mass of the dust grain, 
as it moves through the surrounding gas. At this stage, we 
can make the usual simplifying assumption that for a single 
dust grain, the radiation pressure and the drag force terms 
in Eq. ([g]) dominate completely over the other terms (e.g. 
Lamers fc Cassinelli||1999 1 and balance each other. We can 
then write, 



Tva^QrpL^ 
4nr^cmd 



md 



-Id = 0. 



(10) 



We can now write the radial and azimuthal components of 
Eq. separately. Each of these components must identi- 
cally vanish. Thus we get for the radial component. 



Tva^QrpL^ 



— fo: 



(11) 



where is the radial component of the drag force, defined 



+ ai 



(12) 



where ath is the thermal speed given by ath = ^/ikTlii 
and pm„ is the mean molecular mass of the gas. Typically 
for AGB stars with nearly solar abundance with pulsational 
shocks that extend the density structure we can use p « 1.3 
(e.g. |Bowen|T988t . From Eq. ([lOl we can immediately write 
the azimuthal momentum equation for the dust grains as. 



fV 
J D 

md 



= 0, 



(13) 



where the azimuthal component of the drag force, im- 
plying that there is no drag in the azimuthal direction. That 
is, the dust is co-rotating with the gas. Having obtained ex- 
pressions for the radial and azimuthal components for the 
dust grains we can now re-visit the Euler equation for the 
gas. The radial momentum equation for the gas can be re- 
cast in the form. 



I dp GM-,{1-Td 



du 

dr r ' p dr pc 
wherein, we have the usual definition that 

_ rid TTO^QrpL^ 



p AncGMt 



{JxB)r=Q, (14) 



(15) 



Presently, the azimuthal momentum equation for the gas 
can be written, after some re-arrangement, as. 



dr 



r'^Br-^{rB^) = 0. 



dr 



(16) 



It is to be noted that the azimuthal component of the 



drag force does not appear in Eq. ( 16 1 as it vanishes (see 
Eq. (13l). Moreover, since the dust-to-gas ratio is small, 
i.e., ridmd <C p and since the dust and gas velocities are 
expected to be on the same order with the dust velocity ex- 
ceeding the gas velocity by a reasonable fraction of the gas 
velocity, it is reasonable to make the approximation that 
pur^ 3> ridmdvr^ . This allows us to further make the ap- 
proximation that pur'^ = constant in Eq. (j4]). Therefore, we 
can immediately recover the usual expression for the total 
specific angular momentum of the wind as, 

Br 



-rB^ = L = Q.rA, 



(17) 



A-KpU 

where is the Alfven radius, defined as the distance from 
the centre of the star at which the radial magnetic energy 
density is equal to the kinetic energy density, i.e., \pu'^ ~ 

. At this stage, we assume a polytropic equation of state 
for the gas as given by. 



.^0 



(18) 



where, po and po are the pressure and density, respectively, 
at the surface of the star and 7 is the polytropic index. Then 
writing the density and the Lorentz force term of Eq. ( 14 1 



as functions of radius and the radial velocity, we can easily 
obtain an expression for the gas velocity gradient as. 



dw w N{w, x) 
dx X D{w, x) 



(19) 



where, w = u/ua is the gas speed normalised using the 
Alfven speed and x — r/rA, is the radial distance expressed 
in units of the Alfven radius. The quantities N{'w,x) and 
D{w,x) are the numerator and denominator respectively 
and are given by. 



N{w,x) = (2-(ST{wx'^y"^ - 
x{wx'^ -if + S[ix^{w-1] 
and 

D(w, x) — (w^ — jSt^wx'^Y' 



^1 

X 

(l- 



i-Td-eix- xd))^ 

3wx^ + (wi^ -I- l)w) (20) 



[wx 



(wx ) 



In the above equations, the parameters St 



- SqX X 
_ 2kTA 



(21) 



Sg 



and Sn 



along with 7 uniquely determine the 
locations of the critical points, and hence the morphology of 



the family of solutions of Eq. (191. The critical points are. 
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as usual, defined as the locations at which the both the nu- 
merator and denominator vanish, thereby keeping the right- 



hand side of Eq. ( 19 1 finite. The presence of the Heaviside 



function in Eq. ( 20 1 represents the formation of dust at the 
location x = Xd, the dust condensation radius in units of 



the Alfven radius. The critical wind solution of Eq. ( 19 1 will 



yield the gas velocity profile and thereby enable determina- 
tion of other dependent variables, such as the dust velocity 
profile (to be discussed below), the Mach number as a func- 
tion of distance from the star, the azimuthal velocity of the 
gas, the azimuthal component of the magnetic field, the tem- 
perature profile and the density structure of the gas in the 
envelope of the AGB star. 

As mentioned above, we can determine the dust veloc- 
ity profile after having determined u, the gas velocity as 



a function of radius, with the help of Eqs. (14 1 and (151 



This allows us to express the drift speed as (e.g. Lamers & 



Cassinelli (1999l) 



(V -Uj + (v - U) ath - ' 



)2 




0. 



(22) 



This yields the solution (after employing Eq. (19l), 

V 1/2 



v{r) = u(r) -I- 



V 



(23) 



The dust grain number density nd, can be obtained from 
Eq. ([5|. In the current study we are not solving for the dust 
velocity simultaneously with Eq. (191. However, since the 



dust-to-gas ratio is already a small number and decreases 
monotonically in the wind away from the star to become 
even smaller, we shall therefore make the simplifying as- 
sumption that ndTUdl p ~ (S), the average dust-to-gas ratio 
in the wind. This essentially implies that Ud falls off faster 
than p, in the wind; this is a reasonable assumption since 
the dust velocity is expected to exceed the gas velocity. 

Finally the energy flux per second per steradian can be 



determined by expressing Eq. ( 19 1 as a total derivative. This 
yields. 



F 
pur'^ 



+ 



— + - 
2 7 



7 PA 

- 1 pA 



7-1 



GM^jl-TdQir-rd)) _ B^Br Or 
r 4:Tvp u 



,(24) 



It is immediately evident upon inspecting Eq. ( 24 1 that the 



Heaviside function will present a discontinuity in the flux at 
the dust formation radius, therefore, in order to preserve the 
constancy of energy flux across the dust formation interface 
at r = Tfj, it becomes necessary to subtract a constant term 
to the energy fiux outside the dust formation interface, such 
that, 



F{r) 



F{r)\ 



constant. (25) 

pur^ pur^ 

This constant is essentially the difference between the energy 
fluxes on either side of the dust formation radius (r = r^) 
and is given by const = GM^Fd/rd- Such a constant term 
effectively redeflnes the gravitational potential, without al- 
tering the dynamics; i.e., its derivative vanishes, since it is a 



constant and thus, it does not change the solution topology 



of Eq. (19 1. Therefore we can write, 



Fjr Td) 
pur'^ 




B^Br fir _^ GM*r„ 



47rp u 



Td 



(26) 



Similarly, for r ^ Td we obtain Eq. ( |24[ ), thus we ensure that 
flux is constant across the dust formation interface. 

This completes our derivation of the governing equa- 
tions for the hybrid wind model. We see that once a solu- 



tion of Eq. ( 19 1 is determined, i.e., the gas velocity profile, it 



in turn determines all the other relevant variables including 



the dust grain velocity as given by Eq. ( 23 1 and the energy 
fluxes given by Eqs. ( 24|26 1. In the following section we shall 
describe the numerical treatment breifly and the results are 
presented and discussed in Q 



3 NUMERICAL DETAILS 

The ordinary differential equation (ODE) for the gas ve- 
locity gradient, Eq. (19 1 was integrated as an initial value 
problem for a range of initial conditions in the w ~ x phase 
space. The domain of integration was a;o ^ a; 5C 5, where 



xo represents the stellar surface. Table [3] summarises all the 
physical parameters employed for a typical AGB star. The 
first step in the numerical procedure is the determination 
of the critical points, this is described below. It is to be 
mentioned in this regard, that we chose the value of the 
polytropic exponent to be approximately mid-way between 
unity and the values employes in current 2D-axisymmetric 
MHD codes, « 1.13 (e.g. [Keppens fc Goedbloed]|1999[ ). A 
value of unity represents an insothermal equation of state 
and since the envelopes of AGB stars may be well-mixed 
due to convection, we chose a value slightly more than unity 
to resemble a sort of effective cooling. 



3.1 Determination of critical points 

The critical points are the locations in the w — x phase space 
at which both the numerator {N{w, x)) and the denominator 
{D{w, x)) in Eqs. ( |20[ | and ( |21[ ), identically vanish. Once the 
values of the parameters St, Sq, Sn and 7 are established, 
we then proceed to solve the system of non-linear algebraic 
equations given by. 



N{ws,xs) = 
N{wf,Xf) = Q 

D{ws,Xs) = 
D{wf,xf)^0, 



(27) 
(28) 
(29) 
(30) 



where, Xs represents the distance from the photosphere (in 
units of va) at which the gas velocity is equal to the lo- 
cal sound speed, Ws (in units of ua)- Similarly, the point 
{wf,Xf) represents the location in the phase space at which 
the kinetic energy density of the gas is equal to the local 

total magnetic energy density, i.e., = — stt ' 

called fast point. The root finding is accomplished using a 

Levenberg-Marquardt medium-scale root finding algorithm 
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Table 1. Summary of the different parameters for modelling an AGB star hybrid wind 



Parameter 


Symbol 


Value and/or Comment 


Mass 


M 


~ 5M0 


Radius 


Ro 


~ 5OOi?0 


Mass loss rate 


M 


~ 1.6 X per year 


Surface magnetic field strength 


Bo 


~ 1 G 


Bulk gas velocity (radial) at the surface 


Uo 


~ 2 X lO^^^iiescO (vanishingly small) 


Surface temperature (effective) 


To 


~ 3000K 


Stellar rotation rate 


n 


~ 2 X 10"^" rad/s 


Surface escape velocity 


fesc.O 


6.19 X 10'^ cm/s 


Polytropic exponent 


7 


1.06 



(e.g. More 19781. Typical tolerances employed were about 
10~^^ in order to ascertain the zeros of the system of equa- 
tions ( 27 - 30 1 . This procedure is carried out for parameters 



St, Sg and 5*^ on both sides of the dust formation interface. 
Across the interface the only change is that, 



So a 



(31) 



where Sq represents the value of the parameter outside of 
the dust formation interface, while 5*^ represents the value 
inside the dust formation interface. The remaining two pa- 
rameters St and Sn are unchanged across the interface. 
Presently, we describe the procedure employed for deter- 
mining the location of the radial Alfven point and the dust 
parameter Fd- 



3.2 Determination of the radial Alfven point and 
dust parameter Fd 

We begin with a set of parameters St, Sg and Sn that 
are chosen arbitrarily; however, with the constraint that, for 
the given set of parameters, a critical solution is not phys- 
ically possible. The rationale being that a purely Weber- 
Davis wind is not possible for an AGB star. This will be ex- 
plained in detail later, when the results are discussed. Once 
the above mentioned parameters are chosen, the remaining 
parameters ua and va are continuously varied for different 
values of F^ until we are able to achieve a physical critical 
solution. The chief criterion for the latter being that the so- 
lution is continuous through the radial Alfven point and does 
not have a kink and is required to originate at the base of the 
wind sub-sonic, pass through all three critical points, viz., 
the sonic point, the radial Alfven point and the fast point 
and subsequently leave the star super- Alfvenic. The reader 
at this stage is referred to the subsection on determination 
of the critical solution for details (see S 3.3.2 1. Following an 
initial guess for the parameters {ua, '"a) with Fd fixed to 
a certain value, they are varied with typical step sizes of 
10~® until suitable values are obtained. Once ua and rA are 
determined, the temperature at the radial Alfven point is 
determined according to. 



Ta^To 



UqRq 

UAr\ 



1-7 



(32) 



wherein, the parameters at the base of the wind (subscripted 
with 0) are given in Tabled We now turn our attention to 
the matter of integrating the ODE in Eq. ( 19 1 after having 
determined the above mentioned parameters. 



3.3 Integration of the ODE and determination of 
the critical solution 

3.3.1 ODE Integration 

Integration of the ODE was accomplished with the software 
package ODEPACK employing the subroutine DLSODE us- 
ing backward difference formulae and chord iteration with 



the Jacobian supplied (Hindmarsh 1983 Brown & Hind- 
marsh 1989). Initial conditions were supplied at the begin- 
ning of the integration. Typical error tolerances for conver- 
gence testing that were employed were on the order of 10"^'^ 
for both the absolute and relative errors (see Hindmarsh 



1983 1. For a typical integration over the domain xo ^ x ^ 
a step size of 10~* was employed, resulting in typically 10*- 
l(f function evaluations. A hybrid stellar wind software 
package was specifically developed for the current study 
and this was constructed to be capable of reproducing self- 
consistently the entire family of solutions beginning with an 
arbitrary choice of wind parameters. The entire code takes 
approximately an hour to execute on an AMD Opteron® 
844 1.8 GHz processor. 



3.3.2 The critical solution 

The critical solution is rather unique; it passes through all 
three critical points. In the current study the critical solution 
was determined using a tail procedure. For a given set of pa- 
rameters 7, St, Sg and Sn, forward and backward integra- 
tions were carried out from the points {ws,Xs) and {wj,xj). 
The backward integration from {ws,Xs) was carried out all 
the way to the photosphere of the star and similarly, the for- 
ward integration from {wf,Xf) was carried out all the way 
to the outer boundary of the domain at x = 5. In between 
the backward integration from {wf,Xf) was matched with 
the forward integration from {ws,Xs). An initial guess was 
made for the matching point to be half way between x^ and 
Xf defined by Xm- The tails of the forward and backward 
integrations were terminated at this location and the val- 
ues of the gas velocity w and the velocity gradient ^ were 
compared for the two tails to ensure that the conditions 



(Aw : 



dw 
dx 



dw 



^ 10"^ (33) 
< 10-^ (34) 



were both met. If the conditions in Eqs.(33l and (341 were 



not met, then depending upon the sign of Aw, the matching 
point Xm was appropriately shifted either forward or back- 
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ward by a small amount, typically by Axm = 10~ and the 
tail procedure was re-iterated. Once the tail procedure was 
successful the critical solution was considered to be deter- 
mined. This ensured that the critical solution was continu- 
ous through the radial Alfven point. The tolerance employed 



0.3 



in Eqs. (33 1 and (34 1 was considered to be sufficient given 
the fact that the integration step size was Aa; = 10~*. For 
achieving a higher tolerance, a further reduction in the step 
size was found to be necessary, rendering the the procedure 
needlessly lengthy and increasingly cumbersome in a com- 
putational sense. 

This completes our discussion of the numerical details 
regarding determining the complete family of solutions of 



Eq. (19 1. The results are presented in the following section 



and are discussed therein. 



4 RESULTS AND DISCUSSION 

We present the results obtained by integrating the ODE in 
Eq. (191 to obtain the family of solutions. In the current 



study, the following methodology was employed for calcu- 
lating the hybrid wind. First, for a certain set of arbitrary 
parameters {St, Sg, Sn and 7}, Eqs. ( 27|30 \ were solved 
to obtain the location of the slow and fast points, with the 
radial Alfven point located at (1, 1) in the w — x phase space. 
At this stage, the dust parameter Td was set equal to zero, 
meaning that there hasn't been any dust formation in the 
gas. The parameters are chosen such that a pure WD wind 
is not a physical one, i.e., it is not continuous through the 
radial Alfven point (see Figure [s] and discussion thereof). 
The rationale being that for AGB stars, it is not possible 
to have a mass efllux without dust formation in the enve- 
lope, therefore a pure WD wind is explicitly required to not 
be a physical solution. Second, the dust parameter was set 
equal to a fraction such that ^ ^ 1 and the proce- 
dure described in §3.2| for determing the Alfven velocity and 
radius is carried out in tandem with solving Eqs. ( 27|30 1 
to obtain the sonic point and the fast point. Each time this 
yields a set of parameters {ma, ta, {ws,Xs), (w/,a:/)}. With 
these parameters, Eq. ( 19 1 is integrated to obtain the criti- 



cal solution according to the procedure described in §3.3. 2| 
If success is achieved in finding the critical solution then the 
iterations are ceased and the resulting parameters are fixed 
for the given value of F^j. 

We then determine the temperature profile using the 
velocity profile of the critical solution. This is achieved using 
the prescription. 



T(r) = Ta 



uat'a 



(35) 



where, Ta is the Alfven temperature and is given by 
Eq. (32 I. Once the temperature profile is known, it is pos- 
sible to invert Eq. (35 1 to yield the radius rd, at which the 



temperature falls below the dust condensation temperature 
Td- Once the dust condensation radius (rd) is determined, we 
then proceed to determine the family of solutions to Eq. (191 
such that integration from the photoshpere at r = _Ro to 
r = is carried out with Td ~ and integration from 
r = rd to the outer boundary at r = 5rA is done with 
Fd 7^ 0. The same is then done for the critical solution as 
well. 



0.2 



0.1 




150 



Figure 1. Family of solutions of Eq. l |19| with parameters ua ~ 
0.09tiesc,0i ^A ~ 63.93-Ro, F^ = 0.3 and remaining parameters as 
given in Table 1. The dashed line at r 8.65iio represents the 
dust formation radius. 



0.3 




50 100 
r/R„ 



150 



Figure 2. Illustration showing the locations of the different types 
of solutions of Eq. ( |19| l 



As mentioned before, it is implicitly assumed, for the 
sake of analysing a patently simple model, that beyond the 
dust condensation radius the value of Fd is constant. Concor- 
dantly, all the dust forms at the condensation radius and the 
dust-to-gas ratio is small and given by Eq. ([5| . The family of 
solutions determined using the procedures described above 
are shown in Figure [T] For the solutions shown therein, the 
dust parameter was fixed at Fd = 0.3. It can be seen that 
the dust condensation radius is located at rd « 8.65i?o- The 



temperature at this location was determined using Eq. (351, 
to be approximately 1200K. The dust condensation tem- 
perature was chosen somewhat arbitrarily to lie between 
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1000 — 1500K. Within the dust condensation radius Xd, the 
wind solution to Eq. (191 is a purely Weber-Davis-type of 



wind, and thereafter it is a hybrid wind with dust grains in- 
cluded. The topology of solutions in Figure 1 looks typically 
like a WD-solution, with the three critical points clearly vis- 
ible, viz., the sonic point, the radial Alfven point and the 
fast point. The critical solution emerges from the surface 
of the star sub-sonic, gets accelerated through to the dust 
condensation radius, subsequently passes through the three 
critical points and finally emerges super- Alfvenic at large 
distances. 

At this stage it is convenient to classify the different 
types of solutions in Figure [l] The different types of solu- 
tions and their respective locations in the phase space are 
illustrated qualitatively in Figure [2] Therein, it can be seen 
that the unphysical double- valued solutions to the left of the 
sonic point are referred to, in the current study, as Type I 
solutions. The failed wind solutions, directly below the sonic 
point, are called Type II. The unphysical multi-valued solu- 
tions that make loops, between the radial Alfven point and 
the fast point, are designated as Type Ilia. The unphysical 
double- valued solutions adjacent to the loop solutions are 
called Type lllb. 

The solitary unphysical solution that passes through the 
sonic X-type singularity, is called the Bondi solution, this is 
labelled as B, in Figure [2] While the unphysical double- 
valued solution that intersects the critical solution, at ex- 
actly the two Alfven points, is designated the Alfvemc so- 
lution; labelled A. Similarly, the critical solution is labelled 
as C. 

The unphysical wind solutions that start at the photo- 
phere with super-sonic velocities, just to the right of the 
Bondi solution and subsequently pass through the radial 
Alfven point and get decelerated to sub-sonic velocities at 
large distances from the star, are designated as Type IVa. 
Meanwhile, the unphysical double- valued solutions that pass 
through the radial Alfven point alone, are designated as 
Type IVb. Finally, the unphysical double-valued solutions, 
in the region between the critical solution and the Alfvenic 
solution, immediately to the right of the fast point, are re- 
ferred to as Type V. 

As mentioned earlier, the dust grains condense from 
the gas beyond r^, where the temperature falls below the 
condensation temperature. Beyond rd, we threfore solve the 
hybrid ODE with a fixed value of Fd ■ Thus at the dust con- 
densation radius, the two types of solutions, dust-free and 
hybrid, must match, in terms of velocity of the gas. This is 
illustrated in Figure [S] where we have expressed the coor- 
dinates using a logarithmic scale to facilitate examination. 
In Figure [3] the red solid line represents a hybrid wind so- 
lution for a dust parameter of Fd = 0.3, while the black 
long-dash-dotted line that passes through the radial Alfven 
point with a kink at {uA,rA), is a pure WD solution with 
Fd = 0. As can be clearly seen, the kink at the radial Alfven 
point indicates that this solution is not physical. This in- 
dicates that it is not possible to have a pure WD wind 
for AGB stars; only with dust formation is it possible to 
achieve an outflow. The Bondi type solutions for both the 
hybrid wind (long-dashed red line) and the pure WD solu- 
tion (black long-dash-dotted line), can also be seen to pass 
through the respective sonic points. Similarly the two fast 
points can also be distinguished clearly for the two types of 
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2500 



2000 



1500 



1000 
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Figure 3. Plausible hybrid wind solutions with parameters ua ~ 
O.OQiiescO, ^A ~ 63.93-Ro, F^j = 0.3 and remaining parameters as 
given in Table 1. The red solid line and black long-dash-dotted 
line intersecting at the radial Alfven point, are the critical solu- 
tions of the hybrid wind model and pure WD wind model, re- 
spectively. The green-short-dash-dotted lines are possible Type I 
wind solutions of a pure WD wind that can leave the star via the 
hybrid critical solution after dust condensation at the intersection 
points of the green-short-dash-dotted lines and the red solid line. 
The blue dashed line is the temperature profile calculated for the 
hybrid wind critical solution. It should be interpreted using the 
secondary axis. 



solutions; the Alfvenic solutions pass through them. It is to 
be mentioned at this juncture, that inclusion of an outward 
force in the wind due to radiation pressure has the effect 
of suppressing both the sonic point and the fast point to- 
wards the photosphere; the sonic and the fast points for the 
hybrid wind (in red) can clearly be seen to lie inside their 
respective counterparts of the pure WD wind, in terms of 
distance from the star's photosphere. The green short-dash- 
dotted lines are solutions of Type I, for a pure WD wind. 
The points of intersection of the red solid line represent- 
ing the critical solution of the hybrid model and the green 
short-dash-dotted lines, represent different dust formation 
radii. The temperatures corresponding to these can be in- 
ferred using the blue line, which represents the temperature 
profile, determined using the critical solution of the hybrid 
wind, via Eq. (351. Thus, if a wind solution starts off at the 



base of the wind sub-sonic and travels through the envelope 
of the AGB star according to (say) the fourth green short- 
dash-dotted line from the left, then it can pass through the 
dust formation radius at approximately lOOOK, just ahead 
of about 10-Ro and leave the star via the red solid line, the 
critical solution of the hybrid wind. We therefore see that 
the unphysical Type I solutions that start off dust-free, can 
intersect the hybrid critical solution and leave the star as 
a hybrid critical wind. However, there is a constraint. The 
last Type I solution that can leave the star as a hybrid wind 
intersects the critical solution in red, at the sonic point. Any 
solution of Type I that intersects the red solid line, after it 
has turned and become double-valued, does not represent a 
physical hybrid wind. This is illustrated in the black long- 
dash-dotted Type I solution that turns and then intersects 
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Figure 4. Critical azimuthal gas velocity as a function of ra- 
dius. The vertical dotted lines indicate the locations of the dust 
formation radius and the the three critical points. 



the red solid line, just ahead of the hybrid sonic point in 
Figure 3. Therefore, all the green short-dash-dotted lines 
are allowed possible solutions of different hybrid winds. 

Figure [3] also shows that if a hybrid wind is to be 
achieved, there are two possibile scenarios. First, the dust 
formation radius can lie within or at most upon the sonic 
point of the hybrid wind solution. Second, it can lie beyond 
or upon the fast point of the hybrid solution, or exactly at 
the radial Alfven point. Only the first possibility places the 
dust formation temperature in the acceptable range for typ- 
ical AGB parameters so we will focus mainly on this class 
of solutions (however, c.f. Figure |9] and discussion thereof). 
While the latter scenarios do represent legitimate mathe- 
matical solutions, it is however unlikely that the dust for- 
mation temperature be significantly lower than about lOOOK 
and concomitantly, that the dust formation radius in AGB 
stars, should lie as far out as nearly 70iio- Since it is expected 
that the dust formation must occur within a few stellar radii 
in AGB stars, the only plausible solutions are therefore the 
Type I green short-dash-dotted lines, that intersect the hy- 
brid critical line in red. 

For the hybrid critical solution shown in Figures[l]and[3j 
it is possible to calculate the azimuthal velocity profile of 



the gas by manipulating Eq. (171. This is carried out and is 



shown in Figure [4] As can be seen the gas velocity profile 
rises sharply to a maximum value close to the sonic point and 
then falls off less steeply than the rise; this is consistent with 
the usual WD-picture. The locations of the dust formation 
radius and the three critical points are also indicated therein. 

In Figure [5] we have plotted the energy flux per second 
per steradian leaving the star, in the equatorial plane, for 
the critical solution, alongside the various components of the 
energy flux. The solid line shows the total energy flux in the 
gas as function of radius; it is constant for the critical so- 
lution. It is the sum of all the other lines on the plot. This 
is calculated via Eqs. ( 24|26 1. The short-dash-dotted line is 



Figure 5. Plot showing the energy fluxes calculated using 
Eqs. | |24|26[ l for the critical hybrid wind solution with parameters 
uj^ ^ O.OQDesc.Oi '"A ~ 63.93i?o, = 0.3 and remaining parame- 
ters as given in Table 1. The vertical dotted lines indicate the dust 
formation radius and the the three critical points. The solid hori- 
zontal line represents the total constant energy flux. From the top 
short-dash-dotted curve shows the magneto-rotational energy, the 
long-dash-dotted curve represents the enthalpy, the short-dashed 
curve shows the variation in kinetic energy and finally at the bot- 
tom the long-dashed line represents the gravitational potential 
energy of the gas. 



the magneto-rotational energy flux, which is the sum of the 



last two terms in Eq. ( 24 1 . As the distance from the photo- 



sphere increases, the radial magnetic field falls of as ~ 1/r , 
however the rotational energy is expected to peak close to 
the sonic point, as can be inferred from Figure |4j where 
the azimuthal velocity peaks. The two competing terms col- 
lectively yield a large positive sum closer to the star, and 
the contribution diminishes quite rapidly, beyond approxi- 
mately the dust formation radius, as can be seen in Figure 
5. The long-dashed line represents the gravitational poten- 



tial energy of the wind, the third term in Eq. (241 and it 



is therefore negative. The short-dashed line is the kinetic 
energy of the wind (the first term in Eq. ( |24[ |) which gradu- 
ally increases away from the star as the wind is accelerated. 
On the other hand, the long-dash-dotted line is the enthalpy 
(the second term in Eq. ( |24[ |) which gradually decreases from 
a maximum value at the stellar surface where the density of 
the gas is expected to be the highest and follows a power 
law, of the form ~ p^~^ ■ The decrease is quite gradual as 
the exponent is 0.06 since we used 7 = 1.06. The contribu- 
tion of the kinetic energy and the enthalpy are quite small 
in comparison to the magneto-rotational energy. However, 
when the magneto-rotational energy is added to the grav- 
itational potential energy, then the sum is comparable to 
the kinetic energy. Again, the dust formation radius and the 
three critical points are indicated in the figure. It is also to 
be mentioned that, since it is possible to recover the momen- 



tum terms of Eq. ( 14 1 by differentiating Eq. ( 24 1 , therefore 



the slopes of the different lines shown in Figure 5 would in- 
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Figure 6. Dust and gas velocity profiles for two hybrid winds 
with parameters Uj^ 0.09i)esc,0i ''A ~ 63.93-Roi = 0-3 (in 
black) and with parameters Uj^ 0.07vesc,o, ~ 53.05i?o, = 
0.6 (in red). In the former case, the dust-to-gas ratio was (S) = 
1/1000 and in the latter case it was double this value, i.e. (5) = 
1/500. The remaining wind parameters are given in Table 1. The 
dust velocity profile in each case is determined using Eq. \23\. 



dicate the contributions of these terms. In particular, it can 
be clearly seen that the slopes of the magneto-rotational 
energy and the gravitational energy (which includes radia- 
tion pressure) are the most prominent features in Figure 5. 
This indicates that the the most significant contributions to 
accelerating the wind are indeed due to the Lorentz force 
term plus rotational term in Eq. ( |14[ | and the gravitational 
potential which includes the effect of radiation pressure on 
the dust grains. In this regard, it is evident upon inspect- 
ing Figure 5, that beyond approximately the radial Alfven 
point, the gravitational energy flux becomes flat, indicat- 
ing that the acceleration due to the gravitational potential 
modified by radiation pressure, becomes negligible in the 
wind. However, beyond r = va, the magneto-rotational en- 
ergy still has a small contribution and continues to accel- 
erate the wind. In summary, at small distances, {r <« r^), 
both the magneto-rotational terms and the modified grav- 



itational potential terms of Eq. ( 14 \ have a combined and 



pronounced effect in accelerating the wind; however, at large 
distances (r >~ r^) the effect of the gravitational potential 
becomes negligible while there still persists a small contri- 
bution from the magneto-rotational term, that mildly accel- 
erates the wind outward. 

Figure |6] shows the velocity profiles of the dust and the 
gas plotted as a function of distance from the centre of the 
star. The velocity of the dust grains exceeds that of the gas 
bulk, as is required by the hybrid wind model, in order to 
produce an outward drag force on the gas, to accelerate it 
outward. The velocity of the dust grains is determined ac- 



cording to Eq. ( 23 1, once the velocity profile of the gas is as- 



certained, with an assumed average value for the dust-to-gas 
ratio (S). Figure 6 shows the gas and dust velocity profiles 
for two different sets of model parameters. The solid black 



line represents the gas velocity while the dotted black line 
represents the corresponding dust velocity. This model had 
parameters Td ~ 0.3 and (S) — 1/1000. In order to investi- 
gate the effect of changing the average dust-to-gas ratio, we 
kept all other parameters of the model fixed, in particular, 
the radiation pressure mean efficiency and the stellar lu- 



minosity, were kept constant. Then, according to Eq. (151, 



if the average dust-to-gas ratio is doubled then, then ac- 
cordingly, the dust parameter F^ must also double. Thus, 
for the second model's results shown in Figure 6, we took 
Fd = 0.6 and (S) = 1/500. Thus, the red-dashed line repre- 
sents the gas velocity profile, while the red-long-dash-dotted 
line represents the corresponding dust velocity profile. How- 
ever, changing the dust parameter also changes the locations 
of the the three critical points. With Fd = 0.3 the Alfven 
velocity and Alfven radius were found to be ua ~ 0.09uesc,o, 
va ~ 63.93i?o, while for Fd = 0.6, the corresponding values 
were found to be ua ~ O.OSv^scfi, ta ~ 53.05f?o respec- 
tively. It can clearly be seen that within about 507?o (which 
is approximately the location of the radial Alfven point for 
the second model), the dust in the second model (red-long- 
dash-dotted line) has a steeper rise, indicating a larger ac- 
celeration in the wind. Beyond about 50i?o, the acceleration 
of the wind in the second model {Fd = 0.6) starts to decline 
(see Figure 5 and discussion thereof). However, the wind in 
the first model {Fd = 0.3) at this distance, is still getting 
accelerated, therefore it's velocity increases. Thus, when Fd 
is smaller, acceleration due to radiation pressure continues 
to have an effect, out to larger distances from the star. In 
addition, at a distance of about 50-Ro, the temperature in 
the wind is about ~ SOOif for the model with Fd = 0.3 
and about ~ 550/sr for the model with Fd = 0.6. Thus, the 
bulk of the gas is much cooler in the latter case. As a re- 
sult, it is natural that the velocities in the second model 
are slightly lower than the first. That being said, it is to 
be acknowledged that by increasing the dust parameter by 
increasing the average dust-to-gas ratio, the wind gets ac- 
celerated much faster closer to the star. This is expected; 
however, the terminal velocity in this case is lower, as the 
acceleration due to radiation pressure does not have a pro- 
nounced effect out to large distances (r >~ ta)- 

At this stage we turn our attention to the question of 
changing the temperature at the base of the wind. In AGB 
stars, it is likely that due to density pulsations within the 
star, the temperature and density at the stellar photosphere 
are likely to undergo change (e .g.|Lamers fc Cassinelli|1999[ ). 
It has also been suggested by Soker & Clayton ( 1999 1 that 



magnetic-cool spots probably exist in the region around the 
equator of an AGB star, much like our sun. If such is the 
case, then the temperature at the base of the hybrid wind 
will change locally, and this will have an effect on the stellar 
wind beyond the photosphere in the AGB envelope. In or- 
der to investigate this, we constructed a hybrid model with 
an altogether different base temperature. Following |Soker fc] 



Clayton ( 1999 1, we set the base temperature in a magnetic- 



cool spot to be significantly less than the prescribed aver- 
age 3000K, of the stellar photosphere. Accordingly, in order 
to achieve a hybrid stellar wind, that starts with negligible 
velocity at the stellar surface and gets accelerated to super- 
Alfvenic velocities having passed through the dust formation 
radius and the the three critical points, we treated Fd, ua 
and rA as free parameters. As described in §3.2| we varied 
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Figure 7. Plot showing tiie effect of cfianging the temperature at 
the base of the stellar wind. The black lines are for a hybrid wind 
with parameters To = 3000K, Uj^ ^ 0.09vesc,0y fA ~ 63.93-Ro and 
Fjj = 0.3 while the red lines are for a hybrid wind with parameters 
To = 2000K, UA ^ O.OTuesc.o, rA ^ 52.19iJo and Fa 0.62. Both 
winds have identical remaining wind parameters, as given in Table 
1. The dashed lines represent the locations of the respective dust 
formation radii for the two hybrid winds with T^ 1200K, in 
both cases. Dust formation occurs closer to the stellar surface 
when the base temperature is lowered. 



the values of these parameters until a successful critical so- 
lution was achieved. When the base temperature is changed 
to 2000K, we found that in order to achieve a hybrid stellar 
wind, the three free parameters required to take on values; 
Fa « 0.62, UA ~ Omv^scfi and va ~ 52.19i?o. Figure [7] 
shows how the morphology of the solution changes when 
the base temperature is changed. We have merely shown the 
key solutions, namely, the critical solutions, the Bondi and 
the Alfvenic solutions. The solid black lines are the hybrid 
wind solutions when the base temperature is To — 3000K 
and the dust parameter is r^j = 0.3. The lighter red lines 
are for To = 2000K and Fd ~ 0.62. The dashed lines at the 
far left indicate the locations of the respective dust forma- 
tion radii for the two models; in both cases the dust forma- 
tion temperature was taken to be Td ~ 1200K. It can be 
clearly seen in Figure 7, that when the base temperature 
is decreased to 2000K, all the critical points as well as the 
dust formation radius, are suppressed towards the photo- 
sphere. The formation of dust closer to the stellar surface is 
directly related to the steep drop in the temperature profile 
ahead of the photosphere in the AGB envelope. This find- 



ing is consistent with the results of Soker & Clayton ( 1999 1, 
who found that dust formation occurred closer to the stellar 
surface ahead of magnetic-cool spots on the equator of an 
AGB star. Additionally, the stellar wind critical solution for 
To = 2000K has a lower terminal velocity than its counter- 
part for To — 3000K. This is directly related to the fact that 
the bulk of the gas is much cooler when the base temperature 
is lower. If there exist magnetic-cool spots on the surface of 
an AGB star, then the stellar wind properties, namely the 



wind speed and momentum of the outflow will change, ahead 
of the cool spot in the envelope of the star. This will result 
in asymmetric flows between different parts of the star, in 
addition to the fact, that dust formation will occur closer 
to the star. These will directly lead to MHD instabilities, 
that can grow and become unstable and result in asymmet- 
ric mass loss. However, it is to be noted, that investigating 
such effects is beyond the scope of the current study. In or- 
der to investigate instabilities it will be necessary to carry 
out MHD calculation in at least two, if not three dimen- 
sions. By definition, the current steady-state model cannot 
incorporate dynamic effects such as instabilities. Addition of 
effects of radiative transfer, along with internal chemistry, 
that affects dust grain formation, with relaxation of the as- 
sumption of spherical grains, would be the ultimate goal of 
such an endeavour. It has already been shown by |Woitke| 
& Niccolini ( 2005 I using 2-D hydrodynamic codes with ra- 



diative transfer, that it is possible to capture hydrodynamic 
instabilities. The addition of magneto-rotational effects to 
these models would result in richer gas (and dust) dynamics 
and quite different instabilities in the flow altogether, due 
to the presence of MHD efi'ects. The presence of instabili- 
ties in the outflow are the precursors for asymmetric mass 
loss, which has been theorised to be responsible for white 



dwarf kicks (see|Spruit|T998,. Davis et al. 2008 : Heyl|2007a|b 



2008b|a| |Heyl fc Penrice||20d9[ ). In this context, the current 
work is relevant, as it presents a steady-state solution that 
2- and 3-D MHD-dust-driven wind models can reproduce 
by eliminating the time dependence, i.e., setting the d/dt 
terms to zero. It also presents an additional avenue for the 
formation of instabilities, in AGB stellar winds. 

It is also to be noted that, when the dust parameter was 
increased to values closer to unity, it resulted in the critical 
points being suppressed towards the photosphere; this can 
clearly be seen in Figure [7| Concordantly, as the dust pa- 
rameter is increased, it results in dust condensation closer to 
the stellar surface. In Figure [S] we have plotted the location 
of the three critical points as a function of the dust param- 
eter Ffj. The short-dash-dotted line represents the location 
of the sonic point (va) as Td is varied, the solid line shows 
the variation in the radial Alfven point (ta), while the long- 
dash-dotted line shows the change in the fast point (r/). 
Also plotted therein, is the temperature at the radial Alfven 
point; the short-dashed line. As can be clearly seen, in the 
limit of Fd — >■ 1, the three critical points converge and their 
locations approach the stellar surface. The sharp decline in 
the radial Alfven temperature beyond F^ ~ 0.7 indicates 
that the temperature profile in the AGB envelope falls off in 
an extremely steep manner, making the solutions of the gas 
momentum equation implausible and indeed unphysical. 

To determine the dependence of the critical points on 
the dust parameter, we continuously changed the value of F^j 
with a step size of 10"*, within the limits shown in Figure [t] 
and for each given value of Td, we determined the appro- 
priate set of parameters {ua, ta, {ws,Xs) and (wf^Xf)}, 
that yielded a continuous monotonically increasing critical 
solution through the critical points. Following which, the 
temperature at the radial Alfven point was determined ac- 



cording to Eq. ( 32 1 



In simple isothermal dust-driven wind models, a stel- 
lar outfiow is achieved by setting F^ > 1 (e.g. Lamers & 



Cassinelli (19991). This is done to counteract the force of 
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Figure 8. Plot showing tlie effect of changing the dust parameter 
on the morphology of the family of solutions to Eq. | |19| l. In all 
calculations the common parameters for the models are shown in 
Table 1. The short-dash-dotted line shows variation in the sonic 
point (ra) as a function of F^j, the solid line shows the change in 
the radial Alfven point (r^) with changing F^, while the long- 
dash-dotted line represents the change in the fast point (rf) for 
the same case. The short-dashed line traces the dependence of 
the temperature at the radial Alfven point on F^, this should be 
interpreted using the secondary axis. 



gravity and to drive the gas outward. However, in the cur- 
rent study we found that setting the dust parameter to be 
greater than unity, did not yield a set of critical points; i.e., 
they were found not to exist in the domain Rq ^ r ^ 150-Ro, 
for which Eqs. ( |27|30| ) were satisfied simultaneously. Hybrid 
winds were successfully achieved for < Fd < 1. Indeed, as 
is shown in Figure[7| Fd ~ 0.7 is a reasonable physical upper 
limit, where TA(rd ~ 0.7) « 450K and after which point the 
decline in the temperature proceeds very rapidly. It is to be 
mentioned that for all of the calculations carried out to pro- 
duce Figure 7, the temperature at the base of the wind was 
To — 3000K and all the remaining parameters were identical 
to those given in Table 1. 

Finally, for the sake of completeness, we have shown 
in Figure [9] the plausible hybrid wind solutions, should the 
dust formation radius exist outside the fast point. Again 
the parameters of the hybrid wind are identical to those 
shown in Table 1. It is to be mentioned that it is likely that 



dust formation lies within a few stellar radii (e.g. Lamers 
fc Cassinelli| |199"9t; however, since the parameters of Table 
1 may all be scalable to suit an altogether different type of 
star, we have therefore included Figure [9] to complete the 
scenario of dust forming in the envelope or indeed the outer 
atmosphere of the star. 

In Figure [9] the solid black line represents the critical 
solution of a pure Weber-Davis stellar wind, without any 
dust. While the long-dash-dotted red line represents the hy- 
brid wind critical solution with F^ = 0.3 and ua ~ 0.07wesc,o 
and ta ~ 48.20-Ro, the pure WD wind has the same values 
for Ua and ta- As can be seen in this case, the hybrid wind 
is not a physical possibility as it is not continuous through 
the radial Alfven point. On the other hand the pure WD- 




150 



Figure 9. Plausible hybrid wind solutions with dust formation 
occurring beyond the fast point. The red-long-dash-dotted line 
and the black solid line intersecting at the radial Alfven point, 
are the critical solutions of the hybrid wind model and pure WD 
wind model, respectively. The green-short-dash-dotted lines are 
possible Type V wind solutions of a hybrid wind (F^j = 0.3) that 
can leave the star as a dust laden wind after dust formation occurs 
at the intersections with the solid black WD critical solution. The 
thick green solid line represents a possible hybrid wind solution 
with Fa = 2. 



equatorial wind is continuous. The Bondi and Alfvenic solu- 
tions for the two types of solutions, with and without dust, 
are also plotted in long-dash-dotted-red and black-dashed 
lines respectively. The long-dashed lines in red represent the 
hybrid unphysical solutions of Type V. The intersections 
of the pure WD critical solution in black, with the Type 
V hybrid green-short-dash-dotted lines, represent possible 
locations for the dust formation radius. Thus a wind may 
start off as a pure WD wind at the surface of the star, pass 
through all three critical points and then undergo dust for- 
mation beyond the fast point. At this stage the critical so- 
lution may leave the star by following a hybrid solution rep- 
resented by the short-dash-dotted lines in green, after dust 
condensation. The red-long-dashed lines are not plausible as 
they have either turned and become double- valued (the left 
most long-dashed-red dashed line) or they do not intersect 
the solid black line (the long-dashed-red lines to the right 
of the green-short-dash-dotted lines), at least within the do- 
main indicated. The thick solid green line that intersects 
the solid black line at about r « 95-Ro, is a hybrid wind 
solution of Type V. This solution has the dust parameter 
set to a value greater than unity (Td = 2). Thus it can be 
seen that it is possible to have a hybrid wind with F^j > 1 if 
rd > Tf. In this case, the critical points for the hybrid wind 
parameters are unphysical and do not lie within the domain 
Ro ^ r ^ 1507?o (see earlier discussion relating to Figure[8|. 
However in this case, the critical wind solution has already 
been accelerated through the physically possible pure- WD 
critical points before dust condensation occurs in the wind. 
Thus, the hybrid picture can still work with the dust param- 
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eter greater than unity as long as dust condensation occurs 
beyond the fast point. As mentioned before, Figure [9] does 
not apply to AGB stars (since dust condensation likely oc- 
curs within a few stellar radii), it is included here for the 
sake of completeness and understanding the full nature of 
the hybrid wind solutions. 

This completes our discussion of the results of the study. 
In the following section we present our summary and conclu- 
sions along side a short discussion of directions for further 
investigation. 



a successful hybrid wind. This additionally resulted in sup- 
pressing the dust formation radius as well, towards the pho- 
tosphere of the star, consistent with the findings of |Soker] 



& Clayton (19991. Since the velocity of the wind, ahead of 



the magnetic-cool spot in the AGB envelope, was found to 
be appreciably lesser, than the case when the temperature 
was the average equatorial temperature, it was accordingly 
conjectured, that such an effect, would likely produce asym- 
metric outflows, owing to the formation of MHD instabilities 
in the wind. 



5 CONCLUSION 

We present below a brief summary of the paper and there- 
after a short discussion of possible avenues for further work. 



5.1 Summary 

In the preceding discussion we presented a hybrid wind 
model for AGB stars. The model consists of incorporat- 
ing a dust-driven wind with a Weber-Davis MHD equatorial 
wind. The resulting wind momentum equations yielded ex- 
pressions for the radial and azimuthal velocities of the gas 
and the dust. After eliminating the azimuthal equations, two 
radial equations remained in the model that described the 
velocity profiles of both the gas and the dust. In the model 
described in this paper, we explicitly assumed a steady-state 
for the wind dynamics. 

A WD wind was assumed to begin at the surface of the 
star, one that would eventually fail if not for the formation of 
dust grains at a given radius, which allows the hybrid wind 
to leave the star at super-Alfvenic velocities. The dust for- 
mation was assumed to occur abruptly, at a pre-determined 
radius. All the dust grains were assumed to be perfectly 
spherical with identical size. It was implicitly assumed that 
radiation pressure was purely in the radial direction without 
scattering. The opacity of the grains was implicitly assumed 
to be such that all of the radiation impinging on the grains 
was absorbed and imparted momentum to the grains. The 
resulting drag force was assumed to be purely radial as well. 
The rationale was to develop a simple model to delineate 
the key points of the theory. 

The hybrid wind ODE was subsequently solved using 
finite difference methods, for different values of the dust pa- 
rameter Vd- It was found that, in order to achieve a suc- 
cessful hybrid wind, it was necessary for the dust parameter 
to take values such that < Fd < 1, when dust formation 
occurs within the slow point, i.e., rd < rs. The effect of 
changing the dust parameter revealed that F^ > 1, did not 
yield plausible stellar wind parameters. It was found that 
when Fd — >■ 1, all three critical points converged to the the 
stellar surface. 

Finally, the effect of changing the temperature at the 
base of wind was also investigated. The temperature was 
changed from To = 3000K to To = 2000K to represent a 
magnetic-cold spot on the equator of an AGB star. It was 
found that lowering the temperature not only changed the 
morphology of the family of solutions by suppressing the 
critical points towards the stellar surface, but also required 
a greater value of the dust parameter in order to achieve 



5.2 Avenues for further investigation 

The question of MHD instabilities is an intriguing one. It 
presents a direct route for the onset of asymmetric outflows, 
that have been theorised to cause kicks to the nascent white 
dwarf within an AGB star (sc e|Spruit|19 98: Davis ct al. 2008| 
pHeyl 2007a b, 2008b a, Hcyl fc Penrice|2009 ) . Hydrodynamic 
instabilities have already been captured in 2-D simulations 
of AGB winds (see [Woitke fc Niccolini lf2005' ' Woitke|[2005l 
2006[ [Woitke fc Quirrenbach|2008[ [Woit kc 2008). Therefore 



the next logical step would be to incorporate magnetic fields; 
a complicated step, but one that is necessary in order to get 
a more complete picture of these stars. In order to realise the 
onset of MHD instabilities in the flow, a starting point would 
be a 2-D axisymmetric model with magneto-rotational ef- 
fects coupled with dust formation in the envelope. The re- 
sults of the current paper could be used as a check for the 
steady-state solution of such a model. Such an endeavour 
would undoubtedly yield interesting results and would shed 
new light upon the formation of instabilities in the outflows 
from these stars and answer the question of whether such 
instabilities can lead to appreciably asymmetric mass loss 
and momentum transfer. 

A second avenue, would be to relax one, or indeed sev- 
eral, of the assumptions that were made in deriving the cur- 
rent model. As a first experiment, it may be possible to as- 
sume that there also exists drag in the azimuthal direction. 



This would result in a modification of Eqs. ( 13 1 and ( 16 1 



to include an azimuthal drag term. Concomitantly, the as- 
sumption that the dust-to-gas ratio is small can be relaxed so 
that Eq. Q cannot be approximated. This would ultimately 
result in four coupled ODE's (one each for the radial and az- 
imuthal momenta of both the gas and the dust) that can be 
solved simultaneously to yield the dust and gas velocity pro- 
files in both the radial and azimuthal directions. Third, the 
dust parameter can be assumed to vary with radius rather 
than be kept fixed, this may be the easiest to implement. 
Depending upon the nature of the dependence Fd(r), it will 
change the topology of the solution to Eq. ( |19[ ), since the 
solutions of Eqs. ( 27p0 1 will change appreciably. Fourth, 
the dust grain sizes can be assumed to have a distribution, 
rendering the determination of the drag force more tedious, 
but definitely closer to reality. Finally, there can even be 
assumed to exist a certain degree of scattering, which will 
effectively change the radiation pressure term in Eq. All 
these changes will necessarily make the computation more 
intensive, but will all ultimately yield good dividends and 
further the understanding of AGB winds, complementing 
the current state-of-the-art models. 
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